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Fractoluminescence characterization of the energy dissipated dur- 
ing fast fracture of glass 
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PACS 8i.40.Np - Fatigue, corrosion fatigue, embrittlement, cracking, fracture, and failure 
PACS 83 . 50 . -v - Deformation and flow 

PACS 78.55.Qr - Amorphous materials; glasses and other disordered solids 

Abstract -Fractoluminescence experiments are performed on two kinds of silicate glasses. All the 
light spectra collected during dynamic fracture reveal a black body radiator behaviour, which is 
interpreted as a crack velocity-dependent temperature rise close to the crack tip. Crack velocities 
are estimated to be of the order of 1300 m.s -1 and fracture process zones are shown to extend 
over a few nanometers. 



Although glass is considered as the archetype of brittle 
elastic materials, it was shown in various situations that, 
during fracture, a large part of the stored elastic energy 
is dissipated in permanent deformation of the material. 
It was argued to be mostly dissipated in the formation 
of nanocracks in ultraslow stress corrosion cracking condi- 
tions [Tj[2] . In vacuum, Molecular Dynamics (MD) simu- 
lations [3] of dynamic fracture predict that energy dissipa- 
tion results both from plastic deformation and from bond 
breaking ahead of the crack tip, the two phenomena acting 
at the scale of a few nanometers. The latter predictions are 
in agreement with experiments due to S. Wiederhorn [4j, 
who performed accurate measurements of the energies of 
fracture in an inert environment for different glass compo- 
sitions, and showed that they are approximately ten times 
higher than the typical surface energy values estimated 
by Griffith [5]. However, in order to understand fully the 
nature of the dissipation processes, it is highly desirable 
to measure the size of the Process Zone (PZ) where they 
take place for a given crack velocity. Direct measurements 
being hardly imaginable because of the small dimensions 
involved, and because crack velocities are larger than a 
few hundreds of meters per second, one has to resort to 
indirect characterizations. The scope of the present work 
is to show how fractoluminescence can be used for that 
purpose. 



During rapid fracture, the emission of neutral particules, 
ions, electrons and photons has indeed been observed in a 
wide variety of materials 6-8 . Fractoemission has been 
used to probe both fracture mechanisms and fracture sur- 
face chemistry [9] . Since an early paper by Wick et al. 10 



a growing interest has fostered in oxide glass fractolumi- 
nescence, and the first wavelengths spectra were obtained 
in the eighties 11 . 



These spectra usually exhibit both an energy continuum 
and peaks corresponding to discrete energy values. Chap- 



man and Walton 12 performed dynamic fracture experi- 



ments on different glasses, and showed that the continuum 
corresponds to a black body radiator spectrum. Using a 
model developped by Weichert and Schonert 13 , they 



compared the numerical solutions obtained by these au- 
thors to their own observations in order to evaluate the 
size of the heated propagating zone ahead of the crack 
tip, and found it to be of the order of a few nanometers. 
More recently, Gonzalez and coworkers 14 showed that 



no photons were detected for crack velocities smaller than 
10~ 2 m.s -1 in soda-lime silica. 

The central point of this letter is to provide quantita- 
tive evidence for the above scenario by interpreting the 
black body emission spectrum observed during dynamic 



fracture within the framework of Rice and Levy's 15 and 



Freund's [16] models. We estimate both the crack velocity 
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Fig. 1: Experimental setup: (a) Sketch of the DCDC geometry: 
two symmetric cracks propagate from the central hole, and 
the photon sensor is placed on the path of one of them; (b) 
Photograph of the setup. 



and the size of the PZ, where dissipative processes occur 
and the temperature rises. 

Experimental Setup. — Double cleavage drilled 
compression (DCDC) specimens (Fig. [T]) with imple- 
mented precrack obtained by stress corrosion [17] are used 
to perform experiments. Cracks are propagated at a con- 
stant temperature T = (300 ± 2) K under ambient pres- 
sure. Samples made of pure silica (Corning 7980, USA) 
and float glass (Saint Gobain, France) are parallelepipedic 
(5 x 5 x 25 mm 3 , with 10 /zm tolerance). A hole of radius 
a = (500 ± 10) /im is drilled in their center to trigger the 
formation of two symmetric mode I tensile cracks sketched 
in Fig. [l](a) when the sample is submitted to uniaxial 



18 20. To 



compression at imposed displacement speed 
avoid stress corrosion, the specimen is placed in a bath 
of dodecane oil preventing the water molecules penetra- 
tion at the crack tip. This results in dynamic fracture 
as suggested by Wiederhorn 
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Potential light intru- 
sion has been identified and neutralized. In order to check 
that cracks actually propagate within the midplane of the 
sample, we first monitor experiments by imaging the emit- 
ted photons using a large-aperture lens (numerical aper- 
ture: 0.4, focal length: 12mm) onto the 1024x128 pixel 
array (26.6 mmx3.3 mm) of a cooled (193 K) CCD cam- 
era (Andor, DU-401 BR-DD) with a 1:2.75 magnification 
factor. In accordance with standard Charge-Coupled De- 
vice (CCD) operation, the radiation is continuously inte- 
grated and stored as photo-generated charges in the pixels 
of the CCD, until a read-out operation of the accumulated 
charges takes place. In order to reduce the level of dark 
noise, we programmed a read-out every 5 seconds. This 
period has been optimized in order that the accumulated 
dark noise (~ 4 photo-electrons per pixel) remains smaller 
than the signal emitted by a single crack event, whereas 
the probability that the emission by the short-lived crack 
event occurs during a read-out (~ 50 ms) is negligible. The 
detection spectral range is nominally 400 nm -1150 nm. 
Fig. [5](a) and (b) shows two 1.2 x 9.3 mm 2 images which 
are acquired respectively before and during dynamic frac- 
ture. 

In order to measure the photon spectra, a 300-^m- 
diameter fiber bundle is used to collect the emission locally 



on the propagation path on the sample at d = 1.0±0.2 mm 
from the hole as shown in Fig. [l](a). The output of the 
fiber bundle forms the entrance slit of a spectral disperser 
(Andor, Shamrock 163i) coupled to the above-mentioned 
cooled CCD detector. The spectral resolution is 4 nm in 
wavelength. 




Fig. 2: Imaging-mode measurement during dynamic fracture 
on float glass sample. The top and the bottom images are 
acquired respectively before and during fracture. Two cracks 
emerge symmetrically from the hole, starting by being quite 
straight - as observed in (b) - before the specimen shatters. 



Experimental results. — Fig. [3] shows typical spec- 
tra obtained for silica and float glass. The signal inten- 
sity is always observed to be smaller for float glass, for 
which the maximum consistently takes place at lower en- 
ergy values. On the pure silica spectra, we note, on top 
of a continuous spectrum, the presence of a characteristic 
peaks corresponding to nonbridging oxygen hole center re- 
laxation at 650 nm (~ 1.9-eV) [22]. 
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Fig. 3: Experimental spectra obtained on pure silica (red x) 
and on float glass (blue line) during dynamic fracture. The 
450 nm (~ 2.7-eV) and 650 nm (~ 1.9-eV) bands are enlighted 
in gray. 



In order to check that the continuous part of each ob- 
served spectrum is black body type, and to evaluate the 
corresponding temperature, each curve is normalized by 
the spectrum of a thermal source with known temperature 
(Schott KL 1500 LCD). This procedure permits to get rid 
of the overall spectral sensitivity function of the detection 
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system. According to 23 , the Planck distribution func- 
tion in terms of photons number per second and per unit 
wavelength A can be expressed as B\ (T) oc ef ,/ X T_ 1 where 
(3 = hvi/ks — 0.01439 m.K (h,vi and &b being respec- 
tively the Planck constant, the velocity of light and the 
Boltzman constant). The measured renormalized spectra 
were fitted in the wavelength domain 400-700 nm by the 
ratio kB\(T)/B\(Tref), where k is a constant account- 
ing for the difference in duration and spatial extension 
of the sources. This temperature evaluation method was 
first tested on five known temperature spectra obtained 
by the thermal source ranging between 2650 and 3200 K. 
Note that our detector fails to detect IR thermal photons, 
because of the fast thermal diffusion following strongly 
localized heating, and the subsequent temperature decay: 
the IR emission persists longer than the UV-Vis one. Thus 
the fit concerns only the range 400-740 nm in this case. 
The relative error was found to be less than 1% (cf the in- 
set of Fig. [4j. For our glass samples, the fitting range was 
further reduced to 400-600nm because of the existence of 
a peak at 650nm. Temperature rises at crack tips, both in 
silica and float glass samples, are given in Table [T] 
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Fig. 4: Experimental renormalized spectrum ratio 
kB\(T) I B\(Tref) obtained on pure silica glass during 
dynamic fracture. The crack tip temperature is found close 
to T ~ 5000 K - Inset: Example of three known temperature 
spectra ratio obtained by the thermal source (Schott KL 1500 
LCD). 



Local heating induced by non linear deforma- 
tions in dynamic fracture. — The temperature field 
of a crack tip running at a constant velocity v c was de- 
termined by Rice and Levy [15]. They showed that the 
temperature rise within the PZ around the crack tip is 
equal to: 



AT 




(1) 



with K] the stress intensity factor (SIF). v and E are re- 
spectively Poisson's ratio and Young's modulus, while p, 



k and c are respectively the density, the thermal conduc- 
tivity and the specific heat capacity of the material (see 
Table [2). 

In not too fast fracture, i.e. before complex processes 
(such as branching 



occur (such, as branciimg or extended damage for- 
mation 25 ), the crack velocity v c can be related to 
Ki through Freund's elastic continuum mechanics predic- 
tion 16 : 

= 1 - I ' I (2) 
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Ki 



where vr — 3350 m.s -1 is the Rayleigh wave speed in 
glass. Indeed, Sharon and Fineberg 26 validated this 
prediction in glass, for velocities up to v c = 0A2vr ~ 
1400 m.s -1 . For smaller velocities, invoking Eq. |l]) and 
the Dugdale-Barrenblat expression for the PZ size: 



7T /K, 



(3) 



where a y is the yield stress. Several techniques can be 
used to measure it, nevertheless the value is still de- 
bated. The yield stress ranges approximately between 9 
and 12 GPa in silica glass under UHV conditions, depend- 
ing on the loading 
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and references therein. For 
oxide glasses, Wiederhorn estimated a yield stress close 
to 10 GPa [4| w hich is in good agreement both with MD 



simulations 27 28 , where <r y is predicted to be close to 
9GPa, and with the results obtained by Kurkjian on glass 
fibers [29]. Hence, a y is taken equal to 10 GPa in the fol- 
lowing. We can extract an expression for the crack velocity 

aAT 2 



1AT 2 



(4) 



where a = E 2 pck/2vRd 2 (l — v 2 ) 2 . 

Using Eq. Q, we calculate the crack velocity in each 
case. The error bar Sv c on v c can be deduced from the 
error bar 6 (AT) on AT: Sv c /v c = 2S(AT)/ AT(l+v c /v R ). 
This leads to the results reported in Table [l] 

Table 2: Values at 300 K of parameters used in Eq. Q to 
evaluate the dynamic crack velocity v c . 



Oxide glass type 


Silica 


Float 


Young modulus E (GPa) 


72.7 


72 


Poisson's ratio v 


0.16 


0.23 


Density p (kg.m~ 3 ) 


2201 


2530 


Specific capacity c (J.kg _1 .K _1 ) 


703 


880 


Heat conductivity k (W.m^.K" 1 ) 


1.30 


0.937 


Fracture toughness Kj c (MPa.m 1 / 2 ) 


0.794 


0.749 



Discussion. — We have confirmed Chapman and Wal- 
ton's observations [12] : photons emitted during dynamic 
crack propagation in silicate glasses are mostly of thermal 
origin, indicating a high temperature rise in a small zone 
around the crack tip, where dissipative processes occur. 
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Table 1: AT is the temperature elevation in the PZ and v c is the measured crack velocity corresponding. Ki and R c are 
respectively the stress intensity factor and the size of the PZ. I is the extent of the region where the temperature decreases from 
very high within the PZ to its room value far from it. 





AT 


V c 


Kj 


R c 


£ 


Sample 


(K) 


(m.s -1 ) 


(MPa.m 1 / 2 ) 


(nm) 


(nm) 


Silica 1 


4680±130 


1230±94 


0.99±0.03 


3.9±0.2 


0.68±0.05 


Silica 2 


5480±450 


1480±350 


1.06±0.09 


4.4±0.7 


0.56±0.13 


Silica 3 


4980±210 


1330±160 


1.02±0.04 


4.1±0.3 


0.63±0.07 


Float 1 


5180±840 


1550±730 


1.02±0.16 


4.1±1.3 


0.27±0.13 


Float 2 


4100±530 


1170±410 


0.93±0.12 


3.4±0.9 


0.36±0.12 



For pure silica, however, a peak at 650 nm correspond- 
ing to the relaxation of the nonbridging oxygen hole center 
was observed 
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On the contrary, the 450 nm signa- 
ture of the relaxation luminescence of the oxygen deficient 
center was not observed in our experiments, probably be- 
cause of the elevated OH-content of our material (close to 



1000 ppm) (30j[3T 

However, unlike Chapman and Walton [12], we have 
made no a priori assumption on the value of the applied 
stress intensity factor, the process zone size or the crack 
velocity. Using Rice and Levy's model [15], Freund's elas- 
tic prediction 16 of the crack velocity dependence on the 



fracture energy and a Dugdale-Barrenblatt expression of 
the plastic zone size, we have derived a univocal relation 
between the temperature rise and the crack velocity. For 
pure silica and for float glass, dynamic crack velocities 
are found to be lying between 1172 and 1551m. s _1 how- 
ever, values exceeding the limit velocity ~1400m.s~ 1 are 
extrapolated and may not be trusted (Silica2 and Floatl 
samples in Table [I]). 

Moreover, since the crack propagates in the sample mid- 
plane, we can estimate the value of the stress intensity 
factor Kim f° r a crack length equal to the distance d be- 
tween the hole and the detector from the measured load 
to fracture [20]. In our DCDC geometry, the normalized 
crack length d/a = 2 is slightly lower than the domain 
where the normalized SIF equation was adjusted, it is thus 
reasonable to extrapolate for estimate purpose the stress 
intensity factor values: For the three silica samples, we get 
Kf£ 2 = 0.98, 1.03, 1.00 MPa.m 1 / 2 , and in the float glass 
cases Kf^ at = 1.03,0.92 MPa.m 1 / 2 . The uncertainty on 
Ki m being estimated to 0.05 MPa.m 1 / 2 , the macroscopic 
values are in excellent agreement with the instantaneous 
values of Ki estimated independently from the tempera- 
ture elevation. 

Plastic zone sizes were shown to be ~ 3 to 4 nm in all 



cases, which is a reasonable order of magnitude 32 33 



Note however that Molecular Dynamics simulations per- 
formed on pure silica for a crack velocity v c ~ 300 m.s -1 
predict a slightly larger value, close to ~ 10 nm, with a 
yield stress of 9 GPa [3 28 . This can seem in disagree- 



ment with our observations, since in dynamic fracture the 
plastic zone size is supposed to increase with the applied 



stress intensity factor, and hence with the crack velocity. 
This apparent discrepancy can be due to a different way 
of estimating the PZ size. In Rountree et al's work [3|, the 
whole region where flaws appear is taken into account. In 
our case, we measure the size of the region where the tem- 
perature rise is close to its maximum. However, although 
there is actually a temperature profile within the process 
zone, which could be estimated for metallic materials 34 , 
we can show that for glass, the profile is quite abrupt. As 
a matter of fact, what controls the extent of the region 
where the temperature decreases from it elevated value 
within the PZ to its room value far from it is the ratio 
I — k/pcv c of the thermal diffusion coefficient k/pc to the 
crack velocity v c . This length scale can be shown to be 
equal to ~ 0.6 nm and ^0.3 nm for silica and float glass 
respectively. In both cases, I is much smaller than the PZ 
size R c , which means that we can consider that only the 
PZ is at high temperature while the surrounding material 
is at room temperature. 



We are indebted to B. Pant who started these experi- 
ments and to T. Bernard for his technical help. The au- 
thors also thank J. -P. Bouchaud, C. Fiorini, R.K. Kalia 
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